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Abstract 

This paper is concerned with the dynamical properties of determin- 
istically modeled chemical reaction systems with mass-action kinetics. 
Such models are ubiquitously found in chemistry, population biology, 
and the burgeoning field of systems biology. A basic question, whose 
answer remains largely unknown, is the following: for which network 
structures do trajectories of mass-action systems remain bounded in 
time? In this paper, we conjecture that the result holds when the 
reaction network is weakly reversible, and prove this conjecture in the 
case when the reaction network consists of a single linkage class, or 
connected component. 



1 Introduction 

Building off the work of Fritz Horn, Roy Jackson, and Martin Feinberg 
[8j HI [TOj [121 H31 EE] the mathematical theory termed "Chemical Reaction 
Network Theory" has, over the past 40 years, determined many of the ba- 
sic qualitative properties of chemical reaction networks and, more generally, 
models of population processes. As the exact values of key system param- 
eters, termed rate constants and which we will denote by usually 
difficult to find experimentally and, hence, are oftentimes unknown, the re- 
sults tend to be independent of the values of these parameters. In large part 
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motivated by the Global Attractor Conjecture [6], much of the recent atten- 
tion in this field has focussed on which network structures guarantee that 
trajectories are persistent, in that they can not approach the boundary of 
the positive orthant along a sequence of times [TJ [21 [3J HJ El [151 EEl IUJ- 
In this paper we consider a related question: for which network structures 
do trajectories of mass-action systems necessarily remain bounded in time? 
This question is similar to that of persistence in that both force us to con- 
sider extreme behaviors of the species, and, hence, the monomials of the 
dynamical system. Similar to the well known Persistence Conjecture (see 
below), we will conjecture that all trajectories of weakly reversible systems 
with mass-action kinetics are bounded in time. We will prove this conjecture 
in the case when the reaction network consists of a single linkage class, or 
connected component. The methods used in this paper are similar to those 
introduced in p], where the Global Attracor Conjecture was shown to hold 
in the single linkage class case. 

1.1 Formal statement of the problem 

Two of the most basic questions that can be asked about a mathematical 
model for a chemical, or more generally a population, process are (i) must 
all trajectories be bounded in time and (ii) are trajectories persistent in the 
sense of Definition 11.11 below. 

Definition 1.1. For t > denoting time, let x ) be a trajectory to a 
dynamical system in Mr with initial condition xq. A trajectory (f)(t,xo) with 
state space K> is said to be persistent if 

liminf <j)i(t, Xo) > 0, 

t— >oo 

for all % G {1, . . . ,N}, where (fii(t,x ) denotes the ith component of '</>(£, x ) . A 
dynamical system is said to be persistent if each trajectory with non-negative 
initial condition is persistent. 

Thus, persistence corresponds to a non-extinction requirement. Some 
authors refer to dynamical systems satisfying the above condition as strongly 
persistent [20]. In their work, persistence only requires the weaker condition 
that limsup^^ <f>i(t, xq) > for each i G {1, . . . , iV}. 

The following conjecture of Feinberg (see Remark 6.1.E in [9]) is one of the 
most well known in chemical reaction network theory. It pertains to systems 
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whose reaction networks are weakly reversible, or strongly connected (see 
Section [2]), and is intimately related to the Global Attractor Conjecture [BJ. 

Persistence Conjecture. (Version 1) Any weakly reversible reaction net- 
work with mass- action kinetics is persistent. 

In p], it was pointed out that there are really two natural conjectures 
pertaining to weakly reversible reaction networks with mass-action kinetics, 
and that these should be separated. 

Definition 1.2. For t > denoting time, let 0(t,x o ) be a trajectory to a 
dynamical system in M, N with initial condition xq. A trajectory <j)(t,xo) is 
said to be bounded if 

limsup \4>{t, xq)\ < oo. 

t— >oo 

A dynamical system is said to have bounded trajectories if each trajectory is 
bounded. 

Persistence Conjecture. (Version 2) Any weakly reversible reaction net- 
work with mass-action kinetics and bounded trajectories is persistent. 

Boundedness Conjecture. Any weakly reversible reaction network with 
mass-action kinetics has bounded trajectories. 

Clearly, the latter two conjectures would imply the first, which would 
then imply the well known Global Attractor Conjecture (see [BJ E]). Note 
that none of the conjectures make any assumptions on the choice of rate 
constants, which are the natural parameters found in these systems (see Sec- 
tion [2]). The Boundedness Conjecture stated above is quite similar to the 
Extended Permanence Conjecture found in [7J, which conjectures that all 
"endotactic" systems (which include those that are weakly reversible) are 
permanent. Permanence is an even stronger condition than bounded trajec- 
tories in that all trajectories of a compatibility class (invariant manifold), 
regardless of initial condition, must enter a single compact subset of M> . 
The Extended Permanence Conjecture is proven in [7J in the case when the 
system is two-dimensional. In Section we briefly discuss permanence and 
conclude that weakly reversible, single linkage class systems are permanent 
if there is a 5 > for which liminf^oo <pi(t, xq) > 5 for all xq and all i. That 
is, when the system is, in some sense, uniformly persistent. 

Each of the above mentioned conjectures remains open. In recent years 
there has been a great amount of energy aimed at resolving the Persistence 
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Conjecture, and typically that work has focused on Version 2. This focus 
has been quite natural as much of the motivation for the work stemmed from 
consideration of "complex-balanced" systems, see [H [9] , which are known to 
have bounded trajectories. Relatively little attention has been paid, there- 
fore, to the related Boundedness Conjecture, as formally stated above. We 
will refrain from giving an exhaustive background on the work aimed at re- 
solving the Persistence Conjectures, and instead point the interested reader 
to PQ, where such an introduction, including most of the relevant references 
related to persistence and the Global Attractor Conjecture, can be found. 

1.2 Results in this paper 

In this paper we will prove that the Boundedness Conjecture holds for all 
systems whose network consists of a single linkage class, or connected compo- 
nent (see Section [2]). To prove our results, we will use a method, introduced 
in pQ, for partitioning the relevant monomials of the dynamical system along 
sequences of trajectory points into classes with comparable growths. This 
will allow us to prove that there is a Lyapunov function which decreases 
along all paths when \x(t)\ is sufficiently large. 

We will prove all of our results in a slightly more general setting than 
mass-action kinetics in that we will allow our rate "constants" to actually 
be bounded functions of time. Results pertaining to systems with such a 
generalized mass-action kinetics are useful as these systems arise naturally 
when a system with standard mass-action kinetics is projected onto a subset 
of the species (see Section 3 of pQ). 

The outline of the paper is as follows. In Section [21 we provide a review 
of the requisite definitions and terminology from chemical reaction network 
theory. In Section [21 we give our main results together with their proofs. 

2 Preliminary concepts and definitions 

Most of the following definitions are standard in chemical reaction network 
theory. The interested reader should see [8] or [11] for a more detailed intro- 
duction. 

Reaction networks. An example of a chemical reaction is 2 Si + S2 — > S3, 
where we interpret the above as saying two molecules of type Si combine 
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with a molecule of type S 2 to produce a molecule of type S 3 . For now, 
assume that there are no other reactions under consideration. The Si are 
called chemical species and the linear combinations of the species found at 
either end of the reaction arrow, namely 2Si + S2 and 5*3, are called chemical 
complexes. Assigning the source (or reactant) complex 251 + £2 to the vector 
y = (2, 1, 0) and the product complex S3 to the vector y' = (0, 0, 1), we can 
formally write the reaction as y — > y'. 

In the general setting we denote the number of species by N, and for 
i G {1, . . . ,N} we denote the ith species as Si. We then consider a finite 
set of reactions with the kth denoted by y^ — » y' k , where yk,yi £ ^>o are 
(non-equal) vectors whose components give the coefficients of the source and 
product complexes, respectively. Using a slight abuse of notation, we will 
also refer to the vectors yu and y' k as the complexes. Note that if y^ — 
or y' k = for some k, then the kth reaction represents an input or output, 
respectively, to the system. Note also that any complex may appear as both 
a source complex and a product complex in the system. We will usually, 
though not always (for example, see condition 3 in Definition 12. II below) use 
the prime ' to denote the product complex of a given reaction. 

As an example, suppose that the entire network consists of the two species 
Si and S2 and the two reactions 



where Si — > S2 is arbitrarily labeled as "reaction 1." Then N = 2 and 



Thus, the vector (1,0), or equivalently the complex Si, is both yi, the source 
of the first reaction, and y' 2 , the product of the second. 

For ease of notation, when there is no need for enumeration we will typi- 
cally drop the subscript k from the notation for the complexes and reactions. 

Definition 2.1. Let S = {S,}^ 1; C = {y} with y G Z^ 0; and K = {y ->• y'} 

denote finite sets of species, complexes, and reactions, respectively. The triple 
{S, C, 71} is called a chemical reaction network so long as the following three 
natural requirements are met: 

1. For each Sj G S, there exists at least one complex y G C for which 



Si — > S2 and S2 — > Si, 



(1) 



yi = (l,0), yi = (0,l) 



and 2/2 = (0,1), ^ 2 = (1,0). 
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2. There is no trivial reaction y — > y G TZ for some complex y G C 

3. For any y G C, there must exist a y' G C for which y — )■ y' G TZ or 
y'^yeTZ. 

Notation: We will use each of the following choices of notation to denote 
a complex from C: y, y' , yi, yj, yk, y' k , etc. However, there will be other times 
in which we wish to denote the ith component of a complex. If the complex 
in question has been denoted by yk, then we will write y^. However, if the 
complex has been denoted by y, then we would write its ith component as y i: 
which, through context, should not cause confusion with a choice of complex 
yi. See, for example, condition 1 in Definition 12.11 above. 

Definition 2.2. To each reaction network {S,C,TZ} we assign a unique di- 
rected graph called a reaction diagram constructed in the following manner. 
The nodes of the graph are the complexes, C. A directed edge (y,y') exists if 
and only if y — > y' G TZ. Each connected component of the resulting graph is 
termed a linkage class of the graph. 

For example, the system described in and around ([1]) has reaction diagram 
Si ^ S2, which consists of a single linkage class. 

Definition 2.3. Let {S, C,7Z} denote a chemical reaction network. Denote 
the complexes of the ith linkage class by C C. We say T C C consists of a 
union of linkage classes ifT = U^jLi for some nonempty index set I. 

Definition 2.4. The chemical reaction network {S, C, TZ} is said to be weakly 
reversible if each linkage class of the corresponding reaction diagram is strongly 
connected. A network is said to be reversible if y' — >■ y G TZ whenever 
y^y'e TZ. 

It is easy to see that a chemical reaction network is weakly reversible if 
and only if for each reaction y — > y' G TZ, there exists a sequence of complexes, 
yt,...,y r eC, such that y' y x G TZ, y x y 2 G TZ, ■ ■ ■ , y r _i y r G TZ, and 
y r ^ y G TZ. 

Dynamics. A chemical reaction network gives rise to a dynamical system 
by way of a rate function for each reaction. That is, for each y^ — >• y' k G TZ, or 
simply k G {1, . . . , \TZ\}, we suppose the existence of a function = R yk ^ y ' 
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that determines the rate of that reaction. The functions R k are typically 
referred to as the kinetics of the system and will be denoted by /C, or ]C(t) in 
the non-autonomous case. The dynamics of the system is then given by the 
following coupled set of (typically nonlinear) ordinary differential equations 

x(t) = ^2R k {x(t) ) t)(y' k -y k ) i (2) 

k 

where k enumerates over the reactions and x(t) G 1R> is a vector whose ith 
component represents the concentration of species Si at time t > 0. 

Definition 2.5. A chemical reaction network {S, C, 1Z} together with a choice 
of kinetics K. is called a chemical reaction system and is denoted via the 
quadruple {S,C,1Z,)C}. In the non- autonomous case where the R k can de- 
pend explicitly on t, we will write {S,C,TZ,tC(t)}. We say that a chemical 
reaction system is weakly reversible if its underlying network is. 

Integrating (j2J) with respect to time yields 

x{t) = x(0) + ^2(^J Rk(x(s), s)ds^j (y' h - y k ). 

Therefore, x{t) —x(0) remains within S = sp&n{y' k — y k }ke{i,...,R} 101 a U time. 

Definition 2.6. The stoichiometric subspace of a network is the linear space 
S = spa.n{y' k — yk}ke{i,...,\n\}- The vectors y' k — y k are called the reaction 
vectors. 

Under mild conditions on the rate functions of a system, a trajectory 
x(t) with strictly positive initial condition x(0) G 1R> remains in the strictly 
positive orthant R^ for all time (see, for example, Lemma 2.1 of [18J). Thus, 
the trajectory remains in the relatively open set (x(0) + S) H JR^ , where 
x(0) + S := {z G M N | z = x(0) + v, for some v G S}, for all time. In other 
words, this set is forward-invariant with respect to the dynamics. It is also 
easy to show that under the same mild conditions on R k , (x(0) + S) nl> is 
forward invariant with respect to the dynamics. The sets (a;(0) + S)nM> will 
be referred to as the positive stoichiometric compatibility classes, or simply as 
the positive classes. Note that if each of the sets (x(0) + S) HlR> is bounded, 
then all trajectories of the dynamical system are necessarily bounded also. 
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Therefore, the results of this paper are of interest when each positive class is 
an unbounded set. 

The most common kinetics is that of mass-action kinetics. A chemical 
reaction system is said to have mass-action kinetics if all rate functions R k = 
R yk ^ y ' k take the multiplicative form 

p ( \ „yk,l Vk,2 ^,Vk,N /o\ 

where n k is a positive reaction rate constant and y k is the source complex for 
the reaction. For u G IR> and v G R , we define 

V dcf Vl Vat 

where we have adopted the convention that 0° = 1, and the above is un- 
defined if Ui = when Vi < 0. Mass action kinetics can then be written 
succinctly as R k ( Combining (J2J) and gives the following sys- 

tem of differential equations 

x(t) = Y,K k x(t)y*(y' k -y k ). (4) 

k 

We will generalize the equation @ slightly by allowing each to be a 
bounded function of time. See Definition 2.6 in [7] for a similar definition, 
and [S] for another recent treatment of chemical reaction systems with non- 
autonomous dynamics. 

Definition 2.7. We say that the non- autonomous system {S,C,lZ,)C(t)} 
has bounded mass-action kinetics if there exists an r\ > such that for each 

ke{l,...,\TZ\} 

R k (x,t) = K k (t)x Vk , 

where r\ < Kk(t) < 1/rj for all t > 0. Hence, the vector of concentrations 
satisfies 

x(t) = J2^k(t)x(ty h (y' k -y k ). 

k 

We require some final notation. Let v G M. N for some N > 1, and let 
U C {1, . . . , iV}. We write U[j] for the jth component of U. We then write 
v \u to denote the vector of size \U\ with 

v\u,j = (v\u)j = VU[j] 

for j G {1, . . . , \U\}. Thus, v\u simply denotes the projection of v onto the 
components enumerated by U. For example, if iV = 8 and U = {2,4,7}, 
then for any v G M 8 , v \u — (t> 2 , f 4 , v 7 ) G M 3 . 
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3 Main results 



Recall that for any vectors u, v such that u G M> and v G M we define 

u v = u V i ■ ■ -u V ff , where we use the convention 0° = 1. For completeness, we 
recall the following standard definition. 

Definition 3.1. For any set C, we say {Ti} is a partition of C if each Tj is 
non-empty, [j i Ti = C, and for all i ^ j , Ti HTj = 0. 

The following combination of Definition 13.21 and Lemma 13.31 is a general- 
ization of Definition 4.1 and Lemma 4.2 found in pp. While the generalization 
is not made use of in the current paper, we hope that pointing out that the 
function / below can be nearly arbitrary will be beneficial in future work. 

Definition 3.2. Let C denote a finite set of vectors in R . Let x n G M. N 

denote a sequence of points. For D C M. N with {x n } C D, let f : D x C — > 
M>o. We say that C is partitioned along the sequence {x n } with respect to 
/ if there exists a partition, {Ti}f =1 , ofC, where the Ti are termed tiers, and 
a constant C > 1, such that 

W ifUj, Vk £ Ti for some i G {1, . . . , P}, then for all n 

Qf(Xn,Vj) < f(Xn,Vk) < Cf(x n ,yj). 
(ii) if yj G Ti and yu G Tj +m for some m > 1, then 

fibril Di) 

— — — y oo, as n — >• oo. 

f{x n ,yk) 

When f(x,y) = x y , the case considered in both this paper and J]^, we will 
simply say that C is partitioned along {x n }. 

Note that we have a natural ordering of the tiers: T\ >- T2 >- T3 y ■ ■ ■ y 
Tp, and we say T\ is the "highest" tier, whereas Tp is the "lowest" tier. 

The proof of the following lemma is a slight modification of the proof of 
Lemma 4.3 in [TJ and is omitted. 

Lemma 3.3. Let C denote a finite set of vectors in R Ar . Let x n be a sequence 
of points in R* . For Del" with {x n } C D, let f : D x C -> M >0 . Then, 
there exists a subsequence of {x n } along which C is partitioned with respect 
tof. 
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The following lemma, which is similar in spirit to Farkas' Lemma, states 
that for any set of vectors in M , either their span includes a non-zero vector 
in the non-positive orthant M< , or there is vector normal to their span that 
intersects the strictly positive orthant. 

Lemma 3.4 (Stiemke's Theorem, [19]). Fori = 1, . . . ,n, letui G M m . Either 
there exists a c 6 I" such that 



m 



and such that at least one of the inequalities is strict, or there is a w G M> 
such that w ■ Ui = for each i G {1, . . . ,n}. 

Corollary 3.5. Fori = l,...,n, let Ui G W n . Let U C {l,...,m} and 

V = U c . Either there exists ocGl" such that 




.1=1 

n 

E 

,i=i 



j 



QUi > 0, j G V 



and such that at least one of the inequalities is strict, or there is a w G 
with 

Wj > for j G U 
Wj < for j G V 

such that w ■ Ui = for each % G {1, . . . , n}. 

Proof. Define the vector valued function 9 : M m — > IR m via 



9{u) 



Uj if j G U 
-Uj if j G V 



Applying Lemma [3.41 to the set of vectors 9{uj) proves the result. □ 

Definition 3.6. Let w G M N . The set {i G {1, . . . , N} : w { > 0} is called 
the positive support of w and the set {i G {1, . . . ,N} : Wi < 0} is called 
the negative support ofw. The union of the positive and negative support of 
w, i.e. the set {i G {1, . . . ,N} : Wi ^ 0}, is called the support ofw. 
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Definition 3.7. Let C denote a finite set of vectors in R . Let {T{\ denote 
a partition of C. Let U,V C {1, . . . , N} with U U V nonempty. We will say 
that the vector w e R N is a conservation relation that respects the triple 
(U, V, {Ti}) if the following two conditions hold: 

1. U is the positive support of w and V is the negative support of w. 

2. Whenever yj, yg G T for some i, we have that w ■ (yj — ye) = 0. 

Definition 3.8. Let x n G R> denote a sequence of points. We say that x n 
is partially monotonic if x n ^ > i n+1|i for each i for which lim inf x n ^ = 
and if x n j < x n +i,j f or eac h j for which lim sup^^ x n j = oo. 

Note that Definition 13.81 stands silent on the behavior of those j for which 
< lim inf x nd < lim sup^^ x nJ < oo. 

Theorem 3.9. Let C denote a finite set of vectors in R . Let x n G de- 
note a partially monotonic sequence of points for which lim^oo x n ^ G {0, oo} 
for at least one i G {1, . . . , N}. Let 

U = {ie{l,...,N} : lim^^O} 

n— >oo 

V = {j G {1, . . . , N} : lim x n ■ = oo}. 

n— >oo 

Finally, suppose that C is partitioned along {x n } with tiers Ti, for i = 
1, . . . , P, and constant C > 0. Then, there is a conservation relation w G R^ 
that respects the triple (U, V, {Tj}). 

Proof. We suppose, in order to find a contradiction, that there is no conser- 
vation relation that respects the triple (U, V, {Tj}). Define the sets Wi C R^, 
for % = 1, . . . , P, and W C R N via 

p 

Wi = { Vj -y k eM N \ y„y k G T}, W = \JW i: 

i=i 

and denote the elements of W by {u k }. Note that if T consists of a single 
element, then Wi consists solely of the zero vector. Let m — \U UV\ >0be 
the number of elements in U U V and let Wlt/uv C R m be the restriction of 
W to the components associated with the index set UUV, as discussed at the 
end of Section [21 Denote the elements of V^l^uy by {v k }. Thus, collecting 
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terminology, u k G R , whereas G R m , and for each u k £ W, there is 
a corresponding v k G W^|f/uv f° r which Uk\uuv — v hi however the mapping 
■\uuv need not be injective. 

The set VFlc/uy must contain at least one nonzero vector because other- 
wise any non-negative vector with support U U V would be a non-negative 
conservation relation that respects the triple (U, V, {Tj}), but we have as- 
sumed that no such relation exists. 

Because we have assumed there is no conservation relation that respects 
the triple (U, V, {Tj}), we may conclude by Corollary 13.51 that there exist 
Cfc G R such that 

CkVk - °> if 3 e U 
.v k ew\ uuv I ■ 

v 3 (5) 

^ c k v k > 0, if J G V, 

K v k €W\ uuv ) ■ 

and such that the inequality is strict for at least one j G U U V. 

For v k G W^lc/uv, let m k denote the number of vectors of W that reduced 
to it. Define the function M : R N -> R by 



M{x) 



n 



r x u k \ck/m h 



where and are chosen for u k G if M&li/uy = ^fc G W|i7uv- Note 
that, by construction and by the definition of partitioning along a sequence, 
if u k G W, then there are yj,ye G T, for some i, such that u k = ye — yj and 



1 x w 



(7 — n „fi 



for all n > 1. Therefore, M(x n ) is uniformly, in n, bounded both from above 
and below. Noting that each x n has strictly positive components, we may 
take logarithms and find 



ln(M(x n )) 



(E^V ln ^ ( 6 ) 
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where for a vector u G IR> we define 

ln(w) = (ln(iti), • ■ ■ ,]n(u N )). 
Expanding equation along elements of U U V and (U U V) c yields, 

ln(M(x n )) = Yl 

c k v k \u 1 • ln(x n \u) + I Y c kVk\v ) ■ ln(x n | y ) 



+ ( |^Vuvt) -Hx n \ 



(7) 



By construction, x n ^ is bounded from both above and below for i G (U U V) c . 
Thus, the final term in ([7]) is bounded from above and below. By the inequal- 
ities in ([5]), where at least one term is strict, and the facts that x n> i — > for 
each i G U and x n ,i —> oo for each j G V along this subsequence, we may 
conclude that the sum of the first and second term, and hence ln(M(x n )) 
itself, is unbounded towards positive infinity as n — > oo. This is a contradic- 
tion with the previously found fact that M(x n ) is uniformly bounded above 
and below, and the result is shown. □ 



3.1 Bounded trajectories in the single linkage class 
case 

Define V x : R^ R> by 

N 

K(*)^£>(ln(*) -1) + 1]- (8) 

i=i 

This is the standard Lyapunov function of chemical reaction network theory 
where we have chosen x = (1, . . . , 1) [U [11]. Note that Wi(x) = hxx. It is 
straightforward to show that V\ is convex with a global minimum of zero at 
(!,...,!) j8]. The following is a generalization of Lemma 4.7 in [lj. 



Lemma 3.10. Let {S,C,lZ,JC(t)}, with S = {Si, . . . , Sn}, be a weakly re- 
versible, non- autonomous mass-action system with bounded kinetics. Let 
D C K.> - O ne °f the following two conditions holds: 
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CI: There exists an M > 0, such that for any x G D for which > M or 
Xi < 1/M for at least one i G {1, . . . , N}, we have 



*k{t)x Vk {y' k - y k ) ■ hi(» < 0, for all t > 0. 



k 



C2: There exists a sequence of points x n G D for which lim^oo x n>i G 
{0, oo} for at least one i and 

(i) C is partitioned along x n with tiers {Ti}f =1 , and constant C , and 
(ii) T\ consists of a union of linkage classes. 

Proof. We suppose condition CI does not hold, and will conclude that con- 
dition C2 must then hold. Because condition CI does not hold, there is a 
sequence of points x n G D and times t n > for which lim^oo x n> i G {0, oo} 
for at least one i and 



Applying Lemma 13. 3^ we partition the complexes along an appropriate sub- 
sequence of {x n } with tiers Tj, i — 1, . . . , P, and constant C > 1. Note that 
this also had the effect of only considering the analogous subsequence of {t n }. 
In the following, for tier i G {1, . . . , P}, we denote by 

• all reactions with both source and product complex in T i; 

• {i — > i + m} all reactions with source complex in Tj and product com- 
plex in T i+m for m > 1, 

• {i — > i — m} all reactions with source complex in T, and product com- 
plex in Tj_ m for m > 1. 

Defining u/v = (ui/vi, . . . ,un/vn) for u,v G M> , we re-write the left hand 




(9) 



k 
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side of the inequality OH]) 



^K k (t n )x y n k (y' k - y k ) ■ ln(x n ) 



k 



J2 Kk ^ X n lU 



P r 




) 



(10) 



P-j 



j-1 




(11) 



Note that, by construction, for large enough n any component in the enumer- 
ation ffTT]) is negative, and, in fact, \n.(x V n / x y ^) — > — oo as n — > oo, for these 
terms. The proof that the total summation above (that is, the left hand side 
of (TTTJI) ) must also, for large enough n, be strictly negative unless condition 
C2 holds is now identical to the analogous portion of the proof of Lemma 
4.7 in [TJ, and is omitted here. □ 

Lemma 3.11. Let {S,C,TZ}, with S = {Si, . . . , Sjy}, be a single linkage 
class chemical reaction network. Then, there does not exist a sequence of 
points x n G M> , all in the same stoichiometric compatibility class, for which 
lim^oo x n> i G {0, oo} for at least one i and 

{i) C is partitioned along x n with tiers {Tj}f =l , and constant C , and 
(ii) T\ consists of a union of linkage classes. 

Proof. Note that in the one linkage class case T\ can only consist of a union 
of linkage classes if Ti = C. We suppose, in order to find a contradiction, 
that there is a sequence, {x n }, all in the same stoichiometric compatibility 
class, for which lim^oo x n ^ G {0, oo} for at least one i and 

ii) C is partitioned along x n with tiers {Ti}[ =1 , and constant C, and 

(ii) T\ consists of a union of linkage classes. 

Perhaps after restricting ourselves to a sub-sequence, we may choose x n to 
be partially monotonic (recall Definition 13 .8 j) . Let 



U = {i G {1, N} : lim x 



0} 

oo}. 



V = {j E {1,. . . ,N} : lim x. 



15 



Note that U U V is nonempty by construction. By Theorem 13.91 there is a 
conservation relation w G M> that respects the triple (U, V, {Tj}). 

For each j G V, < and x n j — >■ oo. Thus, if V is nonempty, w ■ x n — >■ 
— oo, as n — ?• oo. If \^ is empty, then {7 is necessarily nonempty and, by 
construction, w ■ x n — > 0, as n — > oo. However, because Ti = C, we have that 
w ■ (y' k — i/k) = for all — > y' k G 7?.. Thus, as the x n are all in the same 
stoichiometric compatibility class, we have that w ■ x n is a constant. This 
shows that we can not have w ■ x n — > — oo, as n — > oo, and so V must be 
empty. However, by our construction we may then conclude that W{ > for 
all i, and Wj > for at least one i. Hence, w ■ x n > 0, and not zero. □ 

We now have our main result. 

Theorem 3.12. Let {S,C,lZ,K,(t)}, with S = {Si, . . . , Sj^}, be a single 
linkage class, weakly reversible, non- autonomous mass-action system with 
bounded kinetics. Then, limsup^^ \<f>(t, Xq)\ < oo for each x ^>o- That 
is, the system has bounded trajectories. 

Proof. Letting D be a non-empty positive stoichiometric compatibility class, 
in the statement of Lemma 13.101 we conclude by combining Lemmas 13.101 
and 13.111 that there is an M > so that for any x G D with |x| > M, we 
have 

K k (t)x Vk (y' k - y k ) ■ ln(x) < 0, for all t > 0. 

k 

Therefore, 

^Vl(0(t,x o )) <0, (12) 

whenever |0(i,x o )| > M. Let B xo = sup{Vx(x) : \x\ = M or x = x }. 
Inequality (|12p shows that Vi(<p(t, Xq)) < B X(i for all t > 0, which when 
combined with the fact that V\(x) — > oo as \x\ — > oo, proves the result. □ 

3.2 Permanence 

Note that Lemma 13.101 and, in particular, equation ( I12p do not give a rate 
at which V\ is decreasing. Hence, we can not conclude in general that all 
trajectories contained within a given stoichiometric compatibility class enter 
a single compact subset of M> . That is, we can not conclude that trajectories 
are permanent in the sense of Definition 13. 131 below. This is quantified above 
by the explicit dependence of B xo upon x . However, we may strengthen our 
results slightly. 
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Definition 3.13. Fort > denoting time, let (f)(t,Xo) be a trajectory to a 
dynamical system in M> N with initial condition Xq. The system is said to be 
permanent if there is a p > such that for every xq G M> , 

p < liminf xo) < limsup <fii(t, xo) < 1/p 
for allie {1,...,N}. 

Lemma 3.14. Let {S,C,lZ,fC{t)} , with S = {Si, . . . , Sjy}, be a weakly re- 
versible, non- autonomous mass-action system with bounded kinetics. Let 
D C JR> be such that dist(D, <9M>o) > 5, for some 5 > 0. One of the 
following two conditions holds: 

CI: For any e > ; there exists an M = > such that for any x G D 
with \x\> M , we have 

K k {t)x Vk {y' k - y k ) ■ ln(x) < -e, for all t > 0. 

k 

C2: There exists a sequence of points x n G D that satisfies lim n _ i>00 \x n \ = oo 
and 

(i) C is partitioned along x n with tiers {Ti}f =l , and constant C , and 
(ii) T\ consists of a union of linkage classes. 

Proof. The proof is essentially the same as for Lemma \'S. 101 We first suppose 
condition CI does not hold. Let e > 0. By our assumption, there must be a 
sequence of points x n G D and times t n > such that lim n _ 5 . 00 \x n \ = oo and 

^2 Kk(tn)x y n(y' k - Vk) ■ ln(x„) > -6. 

k 

The proof is now exactly the same as for Lemma 13.101 except that you 
recognize that for any y G Xi, we necessarily have that x y n — > oo, as n — > oo. 
Therefore, the terms in the summation in (TlT]) not only dominate the others, 
but force the expression to — oo asn-y oo, thereby concluding the proof. □ 

Corollary 3.15. Let {S,C,TZ,K.(t)}, with S = {Si, . . . , S^}, be a single 
linkage class, weakly reversible, non- autonomous mass-action system with 
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bounded kinetics. Let P be a positive stoichiometric compatibility class and 
suppose there is a 5 > so that 

liminf (f>i(t, x ) > 5, for all i G {1, . . . , N} and all x G P. 

t— Kx> 

Then, there is a p > such that for any xq G P 

p < liminf 0j(t, xq) < limsup <fii(t, x ) < 1/p 

for all i G {1, . . . , N}. That is, the system is permanent. 

Proof. The lower bound follows by our assumption. The upper bound follows 
from Lemmas 13.141 (with D equal to P restricted to those x a distance of at 
least 5 away from the boundary). 13.1 1[ and similar arguments as in the proof 
of Theorem 13.121 The only real difference in the proof is that the analog of 
equation ffT2l is 

d 

—V 1 {4>{t,x )) < -e, 

whenever \<f>(t, Xo)\ > M, giving us the needed force to guarantee \(f>(t, Xo)\ 
decreases below some 1/p. □ 

Note that the M = M e> $ > of Lemma 13.141 and hence in the proof of 
Corollary 13. 15[ explicitly depends upon 5. Therefore, it is not sufficient in the 
statement of Corollary 13. 15l to assume the existence of a different 5 = 5 X0 > 
for each xq. Thus, the main results of [1] pertaining to weakly reversible 
networks (with arbitrary deficiency) are not strong enough to guarantee per- 
manence using the above methods. 
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